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LONG-TERM  GOALS 

The  overall  goal  is  to  identify  and  understand  the  physical  processes  that  shape  and  change  coastal 
environments.  Emphasis  is  on  the  application  of  remotely  sensed  signals  that  can  be  compared  with  in 
situ  observations  and  assimilated  within  predictive  models.  In  tidal  flat  environments,  the  major  goals 
are  to  detect  geotechnical  properties  (e.g.,  sediment  strength),  morphologic  features  (e.g.,  channels), 
and  related  hydrodynamic  events  (e.g.,  plumes). 

OBJECTIVES 

The  primary  objective  of  these  joint  efforts  is  to  develop  thermal  methods  for  improved  monitoring 
and  prediction  of  tidal  flat  environments.  Specific  objectives  are  to: 

•  Participate  in  planning  efforts  for  the  Tidal  Flats  DRI,  including  site  evaluation. 

•  Develop  an  integrated  system  for  in  situ  and  remote  (infrared)  measurements  of  thermal 
signals  in  the  field. 

•  Test  and  apply  the  Lovell  [1985]  hypothesis  for  the  porosity  of  sediment  as  a  function  of 
thermal  conductivity. 

•  Explore  inverse  methods  to  optimize  the  assimilation  of  remote  and  in  situ  observations. 

APPROACH 

The  technical  approach  is  to  conduct  laboratory  and  field  experiments  using  simultaneous  remote  and 
in  situ  observations  of  thermal  signals  in  tidal  flat  environments  (Figure  1).  The  experiments  are 
designed  to  study  geotechnical  and  hydrodynamic  aspects  of  tidal  fiats  (Thomson,  N000140710768). 
In  addition,  both  investigators  are  participating  within  the  planning  phase  of  the  Tidal  Flats  DRI 
(Thomson  &  Chickadel,  NOOO 1407 10682),  in  preparation  for  a  series  of  large  collaborative  field 
experiments.  The  planning  phase  includes  meetings  and  field  trips  to  characterize  potential  sites,  as 
well  as  logistical  concerns. 
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Figure  1.  Schematic  diagram  showing  infrared  and  in  situ  measurements  of  thermal  signals  in  a 
tidal  flat  environment.  The  infrared  measurements  of  surface  temperature  are  made  from  an 
elevated  platform,  and  the  in  situ  measurements  of  interior  (both  water  and  sediment)  temperature 

are  made  from  anchored  platforms. 


The  laboratory  and  field  data  will  be  used  to  test  Lovell’s  [1985]  empirieal  formula  for  the  fractional 
porosity  n  (i.e.,  the  water  content)  of  saturated  sediments  as  a  function  of  thermal  conductivity  k,  where 

and  ks,  refer  to  the  thermal  conductivities  of  the  solid  and  fluid,  respectively.  Assuming  a  ID  heat 
balance,  the  temperature  T  at  the  surface  of  the  sediment  (measured  using  infrared  imagery,  see  Figure 
1)  diffuses  downward  in  a  vertical  z  profde  (measured  using  buried  loggers)  at  a  time  t  rate  governed 
by 

dV/dz^  =  {cpik)  dr/dt, 

where  k  is  the  thermal  conductivity  of  interest,  c  is  the  specific  heat,  and  p  is  the  density 
[Subramaniam  and  Frisk,  1992;  Jackson  and  Richardson,  2002],  Sediment  porosity  n  will  be 
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estimated  by  finding  the  best-fit  k  at  eaeh  location  in  the  imagery  and  then  will  be  compared  with 
sediment  samples  collected  by  C.  Nittrouer  &  A.  Ogston  (University  of  Washington). 

The  field  data  also  will  be  used  to  quantify  surface  fluid  velocities  (using  imagery  [Holland  et  al, 
2001])  and  estimate  volume  transport  (using  in  situ  data  [Wunch,  1996]).  These  hydrodynamic 
quantities  will  be  used  to  evaluate  correlations  with  bathymetric  features,  such  as  channels,  and  will  be 
compared  with  velocity  measurements  by  S.  Elgar  &  B.  Raubenheimer  (Woods  Hole  Oceanographic 
Institution). 

WORK  COMPLETED 

During  this  first  year,  investigators  have  attended  DRI  planning  meetings  in  Hawaii  and  South  Korea, 
scouted  field  sites  Washington  State,  begun  development  of  new  methods,  and  collected  preliminary 
data. 
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Figure  2.  Infrared  image  of  a  tidal  flat  (warm  signal)  and  channel  (cold  signal)  atMarisan  Beach, 
South  Korea.  Such  thermal  signals  are  common  on  tidal  flats. 
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The  preliminary  data,  eolleeted  by  Thomson,  include  infrared  images  from  Korean  sites  (Figure  2),  and 
sediment  temperature  profiles  from  Washington  State  sites  (Figure  3).  In  addition,  laboratory  tests  of 
sediment  temperature  loggers  have  been  completed  and  a  new  infrared  imaging  system  is  under 
development.  The  laboratory  tests  validate  the  preliminary  field  measurements  (Figure  3)  and  indicate 
that  the  rms  error  associated  with  the  logging  system,  including  the  effects  of  sediment  anchors,  is  less 
than  0.1  °C.  Processing  of  the  field  measurements  is  ongoing,  including  the  development  of  new 
software  and  methods  to  obtain  best- fit  conductivity  k  values  from  the  data. 


Figure  3.  Tide  (upper  panel)  and  sediment  temperature  profiles  (lower  panel)  recorded  near 
English  Boom  Park,  Skagit  Bay,  WA.  At  low  tide  the  upper  layer  of  sediment  is  heated  by  solar 
radiation,  and  this  heat  is  conducted  downwards  over  time. 

Other  new  methods,  under  development  by  Chickadel,  include  estimation  of  inter-tidal  bathymetry  by 
determining  the  waterline  within  images  at  each  stage  of  the  tide  and  interpolating  to  form  a  Digital 
Elevation  Model  (DEM).  The  DEM  can  be  calculated  at  regular  intervals  to  track  inter-tidal 
morphology  [Ryu  et  ah,  2002].  Eigure  4  shows  an  example  using  infrared  images  from  the 
COFISTREX  experiment  (July  2006),  as  well  as  a  comparison  with  survey  data  (January  2007). 
Infrared  sensing  is  ideal  for  this  method  estimating  bathymetry,  because  the  waterline  can  be  extracted 
day  and  night  and  because  infrared  is  less  prone  to  errors  from  pooled/stranded  water  left  on  the  tidal 
flats. 
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Figure  4.  Inter-tidal  elevation  map  (left panel)  and  transect  (right  panel)  calculated  by  extracting 
waterlines  in  sequential  infrared  images.  The  transect  (colors)  is  in  agreement  with  survey  data 
(pluses)  along  the  black  path  (y  =  114.6  toll  5.3  Km). 


RESULTS 

Results  during  this  first  year  are  primarily  related  to  the  evaluation  of  domestic  field  sites  for  the  Tidal 
Flats  DRI.  Figure  5  shows  satellite  images  of  the  Skagit  Bay  and  Willapa  Bay,  both  of  which  have 
been  documented  thorough  field  trips  and  collection  of  existing  resources  (i.e.,  data,  publications, 
personal  communications).  This  information  has  been  prepared  for  distribution  the  DRI  community  at 
large,  in  support  of  a  draft  science  planning  and  pilot  experiments.  Additional  results  include  the 
confirmation  of  significant  thermal  signals  (eg.  Figures  2  &  3)  at  each  field  site. 

The  Skagit  Bay  site  is  a  predominately  sandy  tidal  flat  with  active  braided  channels,  although  pockets 
of  mud  suggest  that  fine  particles  (i.e.,  mud  and  silt)  are  present  in  the  system  and  being  transported 
through  the  fiats  to  deeper  regions  [McBride  et  al.,  2006].  The  river  input  is  large,  with  an  average 
peak  flow  of  71,310  cfs  (north  and  south  forks  combined, 

http://waterdata.usgs.gOv/nwis/uv712200500).  At  the  northern  end  of  the  bay,  there  is  a  region  of  relic 
(eroding)  mud,  possibly  owing  to  the  construction  of  a  breakwater.  Along  the  southern  perimeter  of 
the  bay,  an  active  depositional  band  of  mud  exists  within  a  few  hundred  meters  of  the  shoreline.  The 
spring  tidal  range  is  4  m  (http://www.cfdnet.com:8080/locations/3237.html). 

The  Willapa  Bay  site  is  predominantly  muddy  tidal  fiat  with  stable  meandering  channels,  although  a 
few  sandy  areas  are  present  at  the  seaward  edges  [Ferraro  and  Cole,  2007].  River  input  is  small,  with 
an  average  annual  peak  flow  of  14,680  cfs  (Naselle  and  Willapa  Rivers  combined, 
http://waterdata.usgs.gOv/nwis/uv712200500,  http://waterdata.usgs.gov/nwis/uv712013500).  The 
southern  portion  of  the  bay  is  extremely  muddy  and  difficult  to  traverse.  The  spring  tidal  range  is  4  m 
(http://www.cfdnet.com:8080/locations/4603.html). 
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Figure  5.  Tidal  Flats  DRI  domestic  field  sites.  The  Skagit  Bay  (left  image)  is  an  approximately  15  x 
15  km  site  centered  at  48°18'52.43  "N  122°27'33.97"W.  The  Willapa  Bay  (right  image)  is  an 
approximately  10 x  10  km  site  centered  at  46°25'14.67”N  123°57'34.46”W. 

IMPACT/APPLICATIONS 

Improving  techniques  to  remotely  quantify  tidal  flat  properties  will  allow  for  real  time  monitoring  and 
safe  operation  in  these  environments.  In  particular,  remote  porosity  estimation  and  channel  detection 
will  improve  navigation  for  amphibious  landings.  In  addition,  the  development  of  techniques  to 
assimilate  remote  and  in  situ  measurements  will  facilitate  the  testing  of  predictive  models. 

RELATED  PROJECTS 

An  ongoing  DURIP  (PI;  Andrew  Jessup)  to  develop  a  “helikite”  aerial  platform  for  scientific  imaging 
will  improve  the  spatial  coverage  and  duration  of  planned  field  experiments. 

An  ongoing  MURI  (Coherent  Structures  in  Rivers  and  Estuaries  Experiment,  PI;  Andrew  Jessup)  has 
provided  infrared  image  data  for  proof  of  concept  applications  in  the  remote  sensing  of  tidal  flats 
(www.cohstrex.apl.washington.edu). 

This  effort  is  a  contribution  to  the  Tidal  Plats  DRI  (www.tidalfiats.org). 
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